Order in the Heisenberg Pyrochlore: The Magnetic Structure of Gd 2 Ti 2 07. 
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The rare earth pyrochlore material GcbT^Or is considered to be an ideal model frustrated Heisen- 
berg antiferromagnet with additional dipolar interactions. For this system there are several untested 
theoretical predictions of the ground state ordering pattern. Here we establish the magnetic struc- 
ture of isotopically enriched 160 Gd2Ti2O7, using powder neutron diffraction at a temperature of 
50 mK. The magnetic structure at this temperature is a partially ordered, non-collinear antiferro- 
magnetic structure, with propagation vector k = 555. It can be described as a set of "q=0" ordered 
kagome planes separated by zero interstitial moments. This magnetic structure agrees with theory 
only in part, leaving an interesting problem for future research. 

PACS numbers: 75.10.-b, 75.25.+Z, 75.50.Ee 
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Lattices based on triangular motifs with antiferromag- 
netic nearest-neighbor exchange are, in the absence of 
further terms in the Hamiltonian, expected to be geomet- 
rically frustrated as all the interactions between coupled 
spins cannot be simultaneously minimized. This situa- 
tion is at an extreme in magnets based on the pyrochlore 
lattice, a face centred cubic array of corner linked tetra- 
hedra. Some thirty years ago Villain considered the 
problem of the Heisenberg antiferromagnet on the py- 
rochlore lattice and argued that it would be a "coopera- 
tive paramagnet" that remained disordered and fluctuat- 
ing down to zero temperature yj] . This has subsequently 
been confirmed by both numerical and analytical calcu- 
lations §-|. However, one question has remained: do 
the dipolar interactions that must be present in a real 
pyrochlore material play a decisive role in its low temper- 
ature physics? In order to answer this question, interest 
has turned recently to the material Gd2Ti207, perhaps 
the only known pyrochlore material that is sufficiently 
free of anisotropy and quenched disorder to realistically 
approximate a dipolar Heisenberg system. 

Confounding expectations, Gd 2 Ti 2 07 was found to 
have an ordering transition at 0.97 K, a temperature that 
is comparable to the Curie- Weiss temperature, 9.6 K Q. 
To understand this transition, two recent theoretical 
studies are particularly important, both of which devel- 
oped the methodology of Reimers et al. |J. The first, by 
Raju et al. [[|, considered a Heisenberg antiferromagnet 
with dipolar coupling and showed that this does not com- 
pletely lift the ground state degeneracy to second order 
in the expansion of the free energy, but rather stabilizes a 
degenerate set of periodic states with h h h propagation 
vectors. It was therefore suggested that the transition 



observed in Gd2Ti207 might be provoked by thermal or 
quantum fluctuations, or alternatively by extra energy 
terms in the spin Hamiltonian. The second study, by 
Palmer and Chalker 0, in contrast, showed that the 
fourth order term in the free energy expansion would 
select, from the degenerate h h h set, an ordered state 
with propagation vector k = 0. In order to test these 
predictions we have now determined the ordered mag- 
netic structure of Gd2Ti207 by powder neutron diffrac- 
tion techniques. 

Natural Gd is almost opaque to neutrons due to the 
high absorption cross section of 155 Gd and 157 Gd, so for 
neutron studies one must use the isotope 160 Gd. Isotopi- 
cally enriched 160 Gd 2 Ti 2 O7 was prepared from 160 Gd 2 O 3 
and Ti02- Stoichiometric quantities were mixed, ground, 
pressed into a pellet and fired at 1200 - 1400 K for ~ 100 
hours. The grinding, pressing and firing processes were 
repeated several times. Powder neutron diffraction data 
were collected using the POLARIS high intensity diffrac- 
tometer at the ISIS pulsed neutron facility. The sample 
was contained in a vanadium can with 3 He exchange gas 
and mounted in an Oxford Instruments dilution refriger- 
ator. Data were collected for 17 hours above and below 
the ordering transition (T c ~ 1 K) at 5 K, 500 mK and 
50 mK. A refinement of the pattern at 5 K showed the 
sample to be phase pure. Below T c ~ 1 K new reflec- 
tions were observed that could be indexed with a mag- 
netic propagation vector of k = \\\- The continuous 
evolution of these reflections as a function of tempera- 
ture showed the transition to be second order. Given 
this observation, representational analysis can be used 
to facilitate magnetic structure determination [p|-pT|. 

We define Gk as the little group of symmetry elements 
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that leave the propagation vector k invariant. The mag- 
netic representation, Tad, of the Gd site (16c in Fd3m), 
can be decomposed in terms of the irreducible represen- 
tations (IRs) of Gk- These are listed in Table | along with 
their associated basis vectors. The application of Gk to 
the four Gd positions of the tetrahedron of the asymmet- 
ric unit results in two orbits. When applied to 0, 0, it 
generates only the position 0,0,0 — orbit 1. When ap- 
plied to another seed position, say , it generates 
the remaining 16c sites: |, |, \ |, |, | and |, |, \ - 
orbit 2. 

These two orbits can be understood in terms of a de- 
scription of the pyrochlore lattice as a set of two dimen- 
sional kagome lattice sheets (the (111) planes) decorated 
by interstitial spins which serve to link the sheets. These 
form 'up' and 'down' pointing corner-sharing tetrahedra. 
If the lattice is described in this way, then orbit 2 de- 
scribes the three atoms that make up the triangular mo- 
tif of a particular kagome sheet and orbit 1 is the unique 
atom present in the interstitial site. The decomposition 
of the magnetic representation Tad for orbits 1 and 2 are: 



j^orbitl_ 
1 Gd - 

Y-orbit1_ 
1 Gd ~ 



or^ + oiY ; + ir 



(i) 



Olf } + 2lf } + Olf } + Hf } + Olf > + 3lf >. (1) 

Landau theory states that in simple systems, a continu- 
ous magnetic ordering transition should involve only one 
IR becoming critical. However, in this case the same 
IRs are not present in the decomposition of T on the 
different sites. For such a situation the only restriction 
that symmetry affords is the constraint that each site, or 
orbit, orders under a single IR with its own set of ba- 
sis vectors. From the decompositions ([!]), there are two 
magnetic representations for orbit 1 and three represen- 
tations for orbit 2 (see Table |). Permutation of these 
results in six allowed combinations of basis vectors, to be 
tested against the experimental data. 

Magnetic structure factors were calculated using the 
GENLES routine of the GSAS suite Q, while the ori- 
entations of the magnetic moments were controlled and 
refined separately by the simulated annealing-based pro- 
gram SARA/i-Refine for GSAS [H. The X 2 and R wp 
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FIG. 1. The intensity of the magnetic reflection as a func- 
tion of the different structural models (i), (ii), (iii) defined in 
the text, for the spin at position (0, 0, 0). 

values for each refinement are shown in Table [n]. It is 
clear that, among the six possibilities, there are only 
two candidates for the magnetic structure; all the other 
combinations can be immediately discarded. We refer to 
these candidate structures as (i) and (ii), see Table ||. 
While both involve the basis vector ipe for orbit 2, they 
differ in that, for orbit 1, structure (i) involves ipi while 
structure (ii) involves the two basis vectors ip2 and ips. 
The basis vector ifig, which refers to the kagome planes, 
in fact corresponds to the "q=0" spin structure observed 
in ordered kagome antiferromagnets |l4f| . For the inter- 
stitial site, structure (i), with the basis vector ipi, corre- 
sponds to a spin pointing parallel to the [1 1 1] direction 
which is perpendicular to the kagome plane. Structure 
(ii), with basis vectors ip2 and ip3, correspond to a spin ly- 
ing with any orientation in a plane parallel to the kagome 



TABLE I. Non-zero IRs and associated basis vectors ip v for the space group Fd3m with k = 555 calculated using the 
program SARA/i- Representational Analysis, jujjn). The labelling of the propagation vectors and the IRs follows the scheme 
used by Kovalev in his tabulated works 
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sheet. We note that powder neutron diffraction cannot 
distinguish this orientation. Three dimensional visualiza- 
tions of the two possible structures are shown in Fig. [j]. 
In both cases the interstitial spins are ferromagnetically 
aligned within their planes, and alternating planes are 
antifcrromagnetically aligned. 

Both structures (i) and (ii) fit well to the data, but 
there is an important difference at » 3.1 A(see Fig. [l]). 
Here, structure (i) predicts more intensity than structure 
(ii) , while experiment suggests that the reflections at this 
d-spacing (for example I, |, |) are systematically absent. 
Thus neither structures (i) nor (ii) are satisfactory. A 
further fit was made which assigned zero moment to the 
interstitial site, structure (iii). This is plausible, since a 
planar antiferromagnetic structure on the kagome sheets 
gives a zero mean field at the interstitial site. 

Structure (iii) , with zero interstitial moment was found 
to improve the fit to experiment by, most importantly, re- 
moving the peak at 3.1 A (see Fig. [l]). The final refined 
profile is shown in Fig. || and goodness of fit parameters 
are listed in Table 0. We therefore conclude that struc- 
ture (iii), with zero interstitial moment, is fully consistent 
with the experimental data. It should be noted that the 
current diffraction experiment cannot shed any light on 
whether the fourth moment is zero as a result of static 
or fluctuating disorder. 

In high symmetry crystals powder neutron diffraction 
does not give an unambiguous determination of the mag- 
netic structure, because there is always the possibility of 
a multi-k structure. In the present case, multi-k struc- 
tures can be formed by contributions from any of the four 
equivalent arms of the star of k. However, recent inves- 
tigations of Gd2Ti 2 07 by Mossbauer spectroscopy EH] 
suggest that all the magnetic spins lie in planes perpen- 
dicular to the [1 1 1] direction below the ordering tem- 
perature. ESR measurements on the system above T c 
have also indicated that the spins lie preferentially in the 
(111) plane JT(|. These observations rule out a multi-k 
structure as the latter would always have a finite spin 
component parallel to [111]. 

TABLE II. Goodness of fit parameters \ 2 and R wp for each 
of the six alternative combinations of IRs, for the data mea- 
sured on the A and C detector banks of POLARIS. The last 
line indicates the fit that considers ordering only in orbit 2. 
The IRs corresponding to models (i), (ii) and (iii) are indi- 
cated. 



Label Orbit A Bank C Bank 





1 


2 


x 2 


Rwp 


x 2 


Rwp 




r 3 


r 2 


43.44 


0.0876 


11.68 


0.0424 


0) 


r. 


r 4 


3.844 


0.0261 


2.720 


0.0205 




r. 


r 6 


13.17 


0.0483 


5.168 


0.0282 




r 5 


r 2 


40.76 


0.0849 


11.62 


0.0423 


(ii) 


r 5 


r 4 


2.838 


0.0224 


2.494 


0.0196 




r 5 


r 6 


11.66 


0.0454 


4.957 


0.0276 


(iii) 





r 4 


2.977 


0.0229 


2.186 


0.0184 




2.0 3.0 4.0 5.0 6.0 7.0 

D-spacing, A 

1 1 1 1 1 1 1 




I I I I I I I 

1.6 1.8 2.0 2.2 2.4 2.6 2.8 
D-spacing, A 



FIG. 2. The final refined magnetic profiles for data col- 
lected from the A and C detector banks of the POLARIS 
instrument from Gd2Ti207 at 50 mK (nuclear reflections are 
marked by the lower line of ticks). 

Our result that k = \\h agrees with the prediction 
of Raju et al. These authors considered a dipolar 
Heisenberg Hamiltonian with coupling parameters ap- 
propriate to Gd2Ti207. They showed that, to quadratic 
order in the Landau expansion of the free energy, the 
ground state consisted of a degenerate manifold of mag- 
netic structures with h h h propagation vectors. It was 
therefore suggested that thermal or quantum fluctua- 
tions, or additional terms in the Hamiltonian, might op- 
erate to select a single ordered state in the real mate- 
rial, with one particular value of h. Palmer and Chalker 
extended this calculation to the fourth order term, and 
found it to be minimized by setting moments of equal 
magnitudes on each site of the elementary tetrahedron. 
A general state was constructed by combining three basis 
vectors of the type ipe (in our notation). Each of these 
basis vectors has three coplanar and one zero moment, as 
described above, but differ in the choice of (1 1 1) plane. 
The condition of equal moments was found to restrict the 
possibilities to a single ground state with h = 0, with the 
four spins in a (1 0) plane, contrary to our findings. 

There are two possible explanations for the disagree- 
ment between experiment and theory. The first is that 



3 



the real material has additional terms in the spin Hamil- 
tonian, such as further neighbour exchange, that over- 
come the fourth order terms in the free energy expan- 
sion. The influence of such terms was considered in a 
rather general sense by Reimers et at || and it would be 
interesting to develop the latter study in the present con- 
text. From an experimental perspective, the Hamiltonian 
could be investigated by studying relatives of Gd 2 Ti 2 07 
such as Gd2Sn 2 07 and Gd 2 SbGa07, where additional en- 
ergy terms might have a different weighting. The second, 
and perhaps more intriguing explanation, is that quan- 
tum fluctuations, neglected in the classical theory, play 
a role in determining the ordering pattern. However, 
the refined value of the ordered moment of Gd2Ti 2 07 is 
6.73 ± 0.05 /is, which is comparable with the maximum 
expected value of 7.94 /i^. This suggests that quantum 
mechanical fluctuations out of the ground state are not 
very important in this spin S = | system. 

In conclusion, the model dipolar Heisenberg antiferro- 
magnet 160 Gd 2 Ti 2 O7 is partially ordered at 50 mK in 
a pattern that consists of "q=0" kagome planes plus a 
zero interstitial moment. It would be interesting to de- 
termine whether this moment exhibits a freezing tran- 
sition at some finite temperature. The partial order of 
chemically pure Gd 2 Ti 2 07 contrasts with the behaviour 
of Heisenberg antiferromagnets with small amounts of 
quenched chemical disorder, such as CsNiCrFg jl7| and 
Y(Sc)Mn 2 These materials show spin freezing tran- 
sitions in the absence of long range order and spin 
correlations characteristic of the pure Heisenberg sys- 
tem B |l7ljl8| . Thus our results support the idea that 
quenched chemical disorder is a strongly relevant pertur- 
bation in the pyrochlore system that is responsible for the 
glassy ground states often observed jl9j. The theory of 
Raju et al. j(| encompasses the experimental behaviour of 
Gd 2 Ti 2 07 in that it predicts an h h h propagation vector. 
We note that, in symmetry terms, h h h is an equivalence 
line in the Brillouin zone, such that all possible ordering 
wavevectors share the IRs and basis vectors that we have 
enumerated in Table H. The selection of h — i might be 
specific to this material or might be more general; this 
question remains to be answered. The higher order the- 
ory of Palmer and Chalker Q does not agree with our 
result. It remains an interesting challenge, both exper- 
imental and theoretical to understand the origin of this 
disagreement. 
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